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Spatially Resolved Images of Dust Belt(s) Around the 
Planet-hosting Subgiant k, CrB 
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ABSTRACT 

Wc present Herschel ^ spatially resolved images of the debris disc orbiting the 
subgiant k CrB. Not only are these the first resolved images of a debris disc orbiting 
a subgiant, but k CrB is a rare example of an intermediate mass star where a detailed 
study of the structure of the planetary system can be made, including both planets 
and planetesimal belt(s). The only way to discover planets around such stars using 
the radial velocity technique is to observe 'retired' A stars, which are cooler and 
slower rotators compared to their main-sequence counterparts. A planetary companion 
has already been detected orbitin g the subgiant k CrB , with revised parameters of 
msinj = 2.1Mj and Upi = 2.8AU ([Johnson et al■|[2008a^ . We present additional Keck 
I HIRES radial velocity measurements that provide evidence for a second planetary 
companion, alongside Keck II AO imaging that places an upper limit on the mass 
of this companion. Modelling of our Herschel images shows that the dust is broadly 
distributed, but cannot distinguish between a single wide belt (from 20 to 220AU) 
or two narrow dust belts (at around 40 and 165 AU). Given the existence of a second 
planetary companion beyond ^3AU it is possible that the absence of dust within ~ 
20AU is caused by dynamical depletion, although the observations are not inconsistent 
with depiction of these regions by collisional erosion, which occurs at higher rates closer 
to the star. 



1 INTRODUCTION 



Our knowledg(| and understanding of exo-planetary sys- 
tems is growing rapidly. Since the first detection of a 
Kuiper-like, planetesimal belt in 1984 (Vega. lAumann et al.l 
19841) )■ the first planet detection around a pulsar in 1992 



Wolszczan fc Frai]||l992l ') and a close-in Ju piter-mass planet 
aroun d a main-sequence star in 1995 l|Mavor fc Oueloa 
Il99d ). the field has exploded. There are now hundreds of 
systems with planet or debris disc detections. There is a 
great deal to be learnt from the growing number of stars 
where both planets and debris discs have been detected. 

Current planet detection techniques are limited to spe- 
cific regions of the parameter space. For example, radial ve- 
locity observations are limited to the inner regions of plan- 
etary systems, whilst direct imaging is limited to the outer 



^ Herschel in an ESA space observatory with science instru- 
ments provided by European-led Principal Investigator consortia 
and with important participation by NASA 



regions. This means that in order to fully characterise a 
planetary system, it is beneficial to have simultaneous ac- 
cess to data from different detection techniques. Radial ve- 
locity observations of A stars on the main-sequence are pro- 
hibited due to hi gh jitter levels and rotationally broadene d 
absorption lines ( Galland et al.ll2005l : iLagrange et al.ll2009l '). 
however, there are now a growing number of detections of 
planets around 'r etired' A stars , now on the subgiant or 
giant branch (e.g Johnson et al.l |2006| . l2007l : iBowler et al.l 
I2OIOI : ISato et al.l I2OI0I ). These provide some key insights 
into the potential differences between the planetary pop- 
ulation around intermediate mass stars, that otherwise can 
only be probed by di rect imaging of planets around main - 
sequence A s tars (e.g.lMarois et al]|2008l: iKalas et al ir2008l') . 
For example, iBowler et al.l ( 2010l ) and lJohnson et al. I (|2010D 
found an increased incidence of giant planets around stars of 
higher stellar mass, as pred icted by planet formation models 
(jKennedv fc Kenvon|[2008l ). 

There are a growing number of sun-lik e stars with 
both planet and debris disc detections (e.g IWvatt et al.l 
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l2012l : iLestrade etlj] |2012| : iLiseau et al] |201C| ). Such sys- 
tems provide key insights into the structure of exo-planetary 
systems and the interactions between planetesimal belts 
and planets. Resolved debris discs often display a va- 
riety of features that can be associated with the pres- 
ence of planets, amongst others, war ps, spirals, brightness 
asymmetries, clumps a nd offsets (e.g . lAugereau et al] I2OOII : 
iMoerchen et alj I2OIII : IWvatt et a"iTll999l 'l. Gaps between 
multiple planetesimal belts could potentially be cleared 
by unseen planetary companions, whilst planets may com- 
monly sc ulpt the inner or outer edges of planetesimal 
belts fe.g.lSu et al.ll2009l:IChiang et aI.ll2009l : IChurcher et alJ 
I2OIII : iLagrange et alj|2012l '). Despite t he ubiquity of debris 
discs aro und main-sequence A stars (|Wvatt et alJ l2007bl : 
iBooth et al. 201 2) and dire ct imaging of a h andful of distant 
planets jMarois et al.ll2008l : ILagrange et aLll2010l ). the inner 
planetary systems remain poorly constrained due to afore- 
mentioned problems with radial velocity measurements. The 
best way to learn about the inner planetary systems of inter- 
mediate mass stars is therefore to observe 'retired' A stars. 
Very little, however, is known about debris discs around such 
'retired' A stars. Such knowledge could act as a further win- 
dow onto the structure of planetary systems around inter- 
mediate mass stars, critical to furthering our understanding 
of planetary systems in general. 

In this work we present Herschel images of a debris 
disc around the subgiant k Coronae Borealis (ft CrB, HD 
142091, HR 5901, HIP 77655) and resolve excess emission 
in the far-infrared, k CrB is a K-type subgiant near the base 
of the giant branch with a m ass of 1.8Mq at a distance of 
31.1pc ljohnson et aljr2008al fl. K CrB is significantly cooler 
than the average main sequence A star, but not significantly 
more luminous, with a lum inosity of 12.3L0 and age of 2.5 
Gyr (| Johnson et aljr2008al ). Radial velocity monitoring o f 
K CrB using the Lick observatories (| Johnson et aLir2008al ) 
found evidence for a planetary companion. The best fit to 
the radial velocity variations find a msini = 2.1Mj planet 
at 2.8 ± Q.IAU, with an eccentricity of 0.125 ± Q.049E|. We 
present far-infrared Herschel observations of this source 
that find and resolve excess emission, alongside follow-up 
radial velocity measurements that suggest the presence of a 
second companion and direct imaging attempts with Keck 
that constrain the potential orbital parameters of this com- 
panion. 

We start by presenting the observations in 33 followed 
by the basic results determined from these observations in 
ij3] Detailed modelling of the Herschel images is presented 
in iJ31 followed by a discussion of the structure of the k CrB 
planetary system in ^and our conclusions are made in ij6l 



2 OBSERVATIONS 

2.1 Keck Radial Velocity Monitoring 

We monitored n CrB at Lick observatory from 2004 to 2009, 
and at Keck observatory from 2010 until present, to search 
for companions to stars more massive than the Sun. This 



Period 


P 


1300 ± 15 days 


Time of pcriccntre passage 


Tp 


13899 ± 160 JD 


Eccentricity 


e 


0.125 ± 0.049 


Argument of pericentre 


OJ 


83.1 ± 29 deg 


Velocity semi-amplitude 


K 


27.3± 1.3ms-i 


Acceleration 


dv 
dt 


1.51 ± 0.52ms"Vr 



Table 1. The new best-fit orbital parameters for k CrB b derived 
from the continued radial velocity monitoring at the Lick and 
Keck observatories. This fit had a reduced chi-squared value of 
1.8 and 7 free parameters, namely, period, eccentricity, longitude 
of periastron, time of periastron passage, global RV ofi'-set, semi- 
amplitude and acceleration. These are derived using the same 
bootstrap Monte Carl o method, as described in further detail in 
Ijohnson et"aLl ll2008al) . 
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^ Calculated using the stellar models of lOirardi et al. I l l2002^ 
2 Updated from the m sin i = l . SMj, 2 .7AU and e = 0.146±0.08 
values quoted in Ijohnson et ah I (jioosi). 
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Figure 1. Radial velocity monitoring of re CrB over 8.09years, 
showing the new orbital fit for re CrB b and a Doppler accelera- 
tion, that provides evidence for a second companion. 

monitoring found the m sin i = 1.8Mj companion at 2.7AU 
l| Johnson et al.ll2008al ) in 2008. Since then, continued mon- 
itoring of this star, over a total of 8.09 years, has updated 
the orbital parameters for k CrB b (shown in Table [1} and 
found msini = 2.1Mj, a semi-major axis of 2.8 ± O.IAU, 
as well as a Doppler acceleration of 1.51 ± 0.52 ms~^yr~^. 
Such a trend provides good evidence for the presence of a 
second companion, however, further monitoring is required 
before the orbit of this companion can be constrained. The 
radial velocity curve for this target is shown in Fig. [1] 

2.2 Adaptive Optics Observations 

Given evidence for the existence of an additional companion 
in the system with a period at least as long as the observa- 
tional baseline (8.09years), ft CrB was observed as part of the 
TRENDS imaging program - a survey dedicated to follow- 
up high-contrast obse rvations of stars showing long-term 
Doppler accelerations (|Crepp et al.ll2012l ). Using NIRC2 (PI 
Keit h Matthews) and the Ke ck II adaptive optics (AO) sys- 
tem jWizinowich et al] |200G| ) , angular differential imaging 
observations were acquired on May 26, 2010 in an attempt 
to directly image the outer body responsible for accelerating 
the star. 

A total of 90 frames were recorded using the narrow 
camera setting. Each frame consisted of a 30 second integra- 
tion time (60 coadds with 0.5 seconds per coadd), resulting 
in a total on-source integration time of 2700 seconds. The 
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Target 


obsID 


Date 


Instrument 


Duration 


K CrB 


1342234353 


15/12/2011 


PACS 


895s 


ft CrB 


1342234354 


15/12/2011 


PACS 


895s 



0.5 1.0 1.5 2.0 

Separation (arcseconds) 

Figure 2. The upper limits on emission from the region sur- 
rounding K CrB determined by the Keck AO imaging. This is 
converted to upper limits on the companion mass in i|3.2l shown 
in Fig. [5] 



field of view (10" x 10", modulo a bad detector quadrant) 
was allowed to rotate to help discriminate between residual 
scattered starligh t (quasi-static spec kles) and faint candi- 
date companions (|Marois et al.ll2C)06l ). The parallactic angle 
changed by 49.5 degrees during the course of the observa- 
tions, allowing us to achieve a close (150 mas) inner-working 
angle. The airmass ranged from 1.04 to 1.08. 

K CrB is sufficiently bright (H=2.58) that observations 
were acquired with the Hcont narrow-band filter. The 300 
mas diameter coronagraphic spot was used to occult the 
star. Unocculted frames were also obtained to measure the 
relative brightness (contrast) between the primary star and 
any off-axis sources. 

Individual frames were processed using standard tech- 
niques to flat-field the array, identify and replace hot pixel 
values, and align and co-add images. We used the locally- 
optimized combination of images (LOCI) algorithm to im- 
prove t he effective signal-to-no ise ratio of speckle suppressed 
frames (|Lafreniere et al.]|2007l ). We did not detect any can- 
didate companions. 

The data, the 5-0" contrasts as a function of angular 
separation, are shown in Fig. [21 These are later converted to 
upper mass limits on the second companion in i]3.2l 



2.3 The Herschel Observations 

Observations were performed using the Herschel Pho- 
todetector and Array Camera & Spectrometer (PACS, 
iPoglitsch erahl (|2010l )') at 100 and 160/im, as listed in Ta- 
ble [2] These observations were performed in mini scan- map 
mode with two observations being performed with a 40deg 
cross-linking angle. Four repeats were used for each observa- 
tion and with eight scan legs per repeat. The total observing 
time was approximately 30 minutes. 

Data were reduced with the Herschel Interactive Pro- 
cessin g Environment version 7.0 Build 1931 (HIPE, lOttI 
l|2010l )) using version 32 of the PACS calibration. Some data 
from the telescope turn-around phase (when scanning above 
5"/s) were used to minimize the ultimate noise level. Maps 
were then made using the HIPE photProject task to provide 



Table 2. The Herschel observations. 



Instrument 


Wavelength 


Photosphere 


Observed 




^m 


mJy 


mJy 


Spitzer 


24 


766 ± 13 


800.1 ±0.1 ±8.0° 


Spitzer 


70 


83 ±2 


426.2'' ±6.5 ±22.3" 


Herschel 


100 


42.34 ±0.69 


335 ± 16mJy 


Herschel 


160 


16.40 ±0.27 


192 ± lOmJy 



Table 3. The photospheric fluxes for k CrB compared to ob- 
serve d fluxes from both archival Spitzer data jKalas &: GrahamI 
I2OO8I) and the Herschel data presented in this work, including 
systematic uncertainties. " For the Spitzer data the first uncer- 
tainties quoted are photometry uncertainties estimated based on 
the data, whilst the second ones include the overall calibration. 
^ The photometry quoted for MIPS 70/im is based on aperture 
photometry rather than the usual PSF fitting photometry, as the 
source is slightly elongated. 



'drizzle' maps (|Fruchter fc Hoo"3l2002l ) with pixel scales of 1 
and 2 arcsec in the 100 and 160 /^m bands respectively. The 
data were high-pass filtered to mitigate low frequency 1// 
noise, using filtering scales of 66 and 102 arcsec (equivalent 
to a filter radius of 16 and 25 PACS frames) in the 100 and 
160 /xm bands respectively. 

The point-spread function (PSF) of the PACS beam 
includes significant power on large scales (10% beyond 1 
arcmin). Consequently, the filtering performed during the 
data reduction will reduce the fiux density of a source 
by 10 — 20%, due to the filter removing the wings of the 
PSF. For point sources this can be readily accounted for 
using correction factors, determined from comparison of 
bright stars with known fluxes with the PACS aperture 
flux. Correction factors of 1.19 ± 0.05 and 1.12 ± 0.05 at 
100 and 160^m were determined from analysis of the DE- 
BRIS (Disc Emission via a Bias-free Reconnaissance in th e 
Infrared / Submilli metre) survey (e.g. Matthews et al-lboiOl ) . 
DEBRIS targets (|Kennedv etatl l2012ai r This can also be 
applied to resolved sources when the source remains similar 
in scale to the beam Full Width Half Maximum (FWHM). 



3 RESULTS 

3.1 Photometry and basic analysis of the Herschel 
data 

By integrating the emission in an area surrounding the cen- 
tral source and comparing this to predictions for the stellar 
photosphere, calculated in the manner discussed in i)4.1l we 
are able to calculate the level of excess flux. These values 
are quoted in Table [31 where the disc is detected at --^ 17 o 
and IScr at 100 and 160^m, respectively. 

Herschel PACS images at 100 and 160/im are shown in 
Figs- Eland 21 Firstly the observations are compared to those 
expected from a single point source, i.e. unresolved emis- 
sion. The residuals after subtracting a point source scaled 
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Figure 3. The 100/im Herschel PACS observations of k CrB. North is up and East is left. The peak position is consistent with the 
stellar position to within 2" (Herschel's 1 — a pointing accuracy), so the observed position is consistent with the Hipparcos astrometry 
projected to December 2011. The colour scale is in mjy/square arcsecond. The hatched circles show the average PACS beam FWHM of 
6.7". The residuals after subtraction of the PSF are shown in the right panel. The contours are at 3, 6, 12, 24, 48, 100 in units of the 
pixel to pixel variation, given by 2.5 X 10~^ or 2.5 X 10~^mjy/arcsec'^ . These residuals clearly show the detection of extended emission 
over and above that of the star. 
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Figure 4. The same as Fig. |3] except at 160/^m. The average PACS beam FWHM is 11.4" and the contours 3, 6, 12, 24, 36 in units of 
the pixel to pixel sigma, given by 6.4 X 10~^ or 2.6 X 10~^mjy/arcsec^ . 



to roughly the peak emission level are shown in the right 
hand panels of Figs. [3] and 21 The residuals clearly show 
that the emission is extended. To learn further about the 
extended emission, we fitted a 2D Gaussian to both sources. 
We find emission that is elongated along a position angle 
of ~ 145°, where the position angle is measured between 
North and the long axis of the ellipse, and East is positive. 
The peak position is consistent with the stellar position to 
within 2" (l-cr pointing accuracy). The full width half max- 
imum of the emission is 11.2" ± 0.01" x 8.26" ± 0.01" at 
lOO^m and 14.9" ±0.013" x 11.3" ±0.017" at 160/im. Given 



that the PSF is extende d by 6.78" x 6.95" at lOO/im and 
12.1" X 10.7" at 160/im (| Kennedy et al.ll2012bl fl. the emis- 



^ These sizes are slightly larger than quoted in the PACS Ob- 
server's Manual, because our data cannot be recentcred on a scan 
by scan basis (i.e. k CrB is much fainter than the PACS calibra- 
tion targets used for PSF characterisation). The PSF sizes are 
quoted as minor x major axis size (i.e. reversed relative to the 
disk sizes) because the PSF tend s to be elongated in the in-scan 
direction l lKennedv et al.ll2012bl) and the k CrB disc is roughly 
perpendicular to the scan direction. 
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sion is clearly resolved at both wavelengths. The elliptical 
shape of such emission is difficult to reconcile with a spheri- 
cally symmetrical shell, but could result from an azimuthally 
symmetrical disc, viewed along a line of sight inclined to the 
disc plane. In which case, the major to minor axis ratio im- 
plies that the disc would be inclined at an angle of 58° ± 1° 
(lOO/xm) or 48° ± 1° (160^m), from face-on, where the un- 
certainties come straight from fitting the 2D Gaussian. The 
disc would have a deconvolved diameter of ~ 280 AU at 
lOO/xm and 320AU at 160^m. 

3.2 Detection limits on (planetary) companions 

The constant acceleration found by the radial velocity mon- 
itoring of K CrB, can be used to place limits on the orbital 
parameters of the second companion to k CrB. The simplest 
assumption to make is that the planet is on a circular orbit, 
in which case a minimum limit on the companion's mass can 
be calculated by assuming that its gravity is responsible for 
the observed acceleration. This limit depends on the separa- 
tion of the companion from the star, or its semi-major axis 
and is given by: 

msini > 1.22M,/ (1) 

where the semi-major axis is a in AU. 

This limit is plotted in Fig. [S] Our non-detection of 
the companion using high-contrast imaging places an upper 
limit on the possible mass of the companion, as a function 
of its projected separation from the star. Our (5a) sensitiv- 
ity to off-axis sourc es, as a function of angular separation, 
calculated using the lBaraffe et al] (|2003l l evolutionary mod- 
els, is over-plotted in Fig. [5] This was calculated using k 
CrB's measmed parallax (vr = 32.79 ± 0.21) and estimated 
isochronal age (2.24 ± 0.15 Gyr). The largest uncertainty 
here is in the age of the primary star, followed by system- 
atic errors in the thermal evolutionary models. 

These non-detection limits complement the radial veloc- 
ity measurements and place tight constraints on the proper- 
ties of K CrB c. The allowable mass and semi-major axis 
parameter space is shown by the shaded area in Fig. [S] 
For instance, Doppler observations rule out the existence 
of extrasolar planets (m < ISMj) beyond ~ 40 AU. The 
minimum possible mass of k CrB c can be determined by 
considering that the minimum possible orbital period of the 
companion, even with an eccentricity of e = 1, corresponds 
to the time baseline of the observations. In this case, the 
minimum period is therefore 8.09 years, which corresponds 
to a minimum semi-major axis of 7.3AU and thus, mass of 
msini > 0.5Mj. 



4 MODELS OF THE Herschel DATA 

The Herschel images of k CrB resemble images of many 
other debris discs around main-sequence stars. The confined, 
disc-like nature of the source suggests that we are observing 
the coUisional remnants of planetesimals in a Kuiper-like 
belt. Although k CrB has evolved off the main-sequence, 
its properties have not changed significantly from that of 
a main-sequence star, in particular, there is no expectation 
that the rate of stellar mass loss has increased beyond the 
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Figure 5. The limits on the planet/companion mass and sepa- 
ration from the host star, based on the non-detections found in 
the AO imaging (black dashed Hne, see i|2.2l l and the radial ve- 
locity trends found in the Keck RV data (solid line- see i|3.2l l. For 
the latter, the mass limit is the minimum mass, msini and the 
projected separations, the semi-major axis of the planet, assum- 
ing that it is on a circular orbit. In order to explain both sets of 
observations, k CrB must have a companion that lies within the 
shaded area. The dotted line shows the minimum planet mass, 
based on the assumption that the planet's orbit is inclined by 
60° from face-on, in the same manner as the disc. 




10' 10° 10' 10- 10' 

wavelength (^m) 



Figure 6. SED for k Cr B. Photometry is shown as black 
dots or black triangles for upper limits. Disk (i.e. photosphere- 
subtracted) fluxes and upper limits are shown as grey dots and 
open triangles. The stellar spectrum is shown as a blue line and 
the modified black-body disk model as a red line, with the total 
shown as a black line. In the modified black- body model Ap was 
arbitrarily set at 210/im and /3 at 1, following IWvatti j2008l) . 

gentle stellar winds of main-sequence stars[f|. Thus, although 
previous explanations for observations of giant stars with 

^ Significant stellar mass loss that could produce an infra-red 
excess emission only occurs towards the tip of the giant branch 
and on the asymptotic giant branches. Evidence against increased 
stellar mass loss rates for sub-giant stars comes from the lack of a 
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infrared excess ha ve include d spora dic mass ejections and 
interstellar clouds (|Kim et"al.. ,2001 : iZuckerman et al.lll995l : 
|juralll999l ). the evidence in the case of k CrB is strongly 
in favour of emission from a debris disc. The only other 
possible source of emission at such wavelengths could be the 
companion(s), which the level and morphology of the excess 
emission lead us to believe is clearly not the origin in this 
case. 

In this section we consider the Herschel observations in 
their own right, comparing the observations with the emis- 
sion from a model debris disc. The intention of this mod- 
elling is to derive as much information as possible regarding 
the disc structure from the Herschel images. Given the lim- 
itations on the information available from these images, it 
will only really be possible to determine a rough size for the 
disc and place some broad constraints on its orientation. 
In ij5] we discuss these models in the light of the limits on 
planetary companions in this system. 

In order to model the emission from a debris disc, 
firstly, the contribution of the stellar emission must be 
accounted for. Optical and near-infr ared photometry is 



19781: Moshir et all 19931: Hauck & MermilhodI 


1997; 


ferryman & 


ESA"l99y: 'Hog ct al."2000l: ICutri et all 


2002; 


Mermilliod 


1987; Ishihara ct al. 20l3). These data 


were 



used to find the be st fitting stellar model, using the 
PHOENIX Gaia gr id (iBrott fc Hauschildl 120051), via a X 
minim isation, as in iKennedv et all (|2012al lbl): IWvatt et al.l 
(|2012h . This method uses synthetic photometry over known 
bandpasses and has been v alidated against high S /N MIPS 
24^m data for DEBRIS (|Matthews et al.1 I2OI0I ') targets, 
showing that the photospheric fluxes are accurate to a few 
percent for main-sequence, AFG-type, stars. 



4.1 Spectral Energy Distribution (SED) 

The synthetic stellar spectrum is plotted in Fig. [B] We have 
added to this the Herschel PACS d ata, as well as archival 
Spitzer data l|Kalas fc Grahamll2008l ) and data points from 
the IRAS faint source catalogue. There is clear evidence for 
excess emission above the predicted stellar spectrum, as can 
also be seen in Table. [S] 

The only information that can be obtained from just 
the spectral energy distribution (SED) is an estimate of 
the disc temperature. In order to determine this, we make 
the simplest possible assumption; that the dust grains emit 
like black-bodies. We anticipate that the inefhcient emis- 
sion properties of real grains reduce the flux at wavelengths 
longer than Aq by a factor (^-^^ ^ , where we have introduced 
the free parameters Ao and /3 to take this into account. Since, 
here we have detections at only four different wavelengths, 
Ao and /? are very poorly constrained, but nonetheless il- 
lustrative of the reduced emission anticipated at long wave- 
lengths, that could be relevant for future observations, for 
example with ALMA. 

The disc temperature can be determined from Fig. [6] 
using our modified black-body description. We find a tem- 
perature of 60 ± lOK, although /3 and Ao remain uncon- 



need to include any stellar mass loss in evolutio nary models that 
fit observations of globular clusters l llbenlll967l ') 



strained and it is clear that some discrepancy exists between 
the IRAS 60/im and Spitzer MIPS 70/im points. It is not 
possible to fit a modified black-body that agrees with both 
these points to within the uncertainties (which are smaller 
than the data points on the plot). We, therefore, assessed 
whether the Spitzer 70fim point might be contaminated by 
background sources, but deem this to be unlikely as the flux 
varies by less than 4% between apertures of different sizes, 
whilst the uncertainty on each data point is 5%. There is 
more reason to question the IRAS point, as there is a greater 
than 30% variation in flux between the IRAS Point Source 
Catalogue (IRAS-PSC) and IRAS Faint Source Catalogue 
(IRAS-FSC), whilst the quoted uncertainties are 10% in 
IRAS-PSC and 8% in IRAS-FSC. In general, the IRAS-FSC 
is more reliable, and this is the point that we use, however, 
it may be that the uncertainty on this point should be in- 
creased above that quoted. 

If this fit is correct and the dust indeed acts like a black 
body, a temperature of 60 ± lOK would correspond to a ra- 
dius of 75I20AU. The large uncertainty in the disk temper- 
ature, and therefore black-body radius, arises because the 
temperature is degenerate with Ao and p. The inferred disc 
radius is therefore smaller than suggested by the images, 
which could arise due to the presence of small (fim) grains 
that are hotter than blackbodies. However, the temperature 
is sufficiently uncertain that this conclusion is not strong. 



4.2 A model of the Herschel images 

K CrB is resolved in both of the Herschel images. We, 
therefore, attempt to fit the observations using a simple 
model. Firstly, the stellar photosphere is subtracted, such 
that we are left with emission from the debris disc, which 
can directly be compared with the emission from a model 
disc. The models we g e nerate here use the m e thods de- 
scribed in IWvatt et all (|2012l ): [Kennedy et~al] (|2012al lbl). 
The basic idea is to determine the spatial distribution of 
the dust, which coupled with the grain emission properties 
at each wavelength can be used to produce a high resolu- 
tion model image of the disc at each wavelength. In these 
models we make no attempts to constrain the grain prop- 
erties or size distribution, since in the absence of spectral 
features this will be degenerate with other assumptions. In- 
stead we use the modified black-body prescription, outlined 
in ^4.11 In this modified black-body prescription, the tem- 
perature of the disc is assumed to be T = /tT);!,, where 
Ttt = 278.3 (L./Lq)* (r/A(7)"i = 521{r/AU)-i is the 
temperature that a black-body at a distance r from the star 
of luminosity L« = 12.3Lq would achieve, and /t is a pa- 
rameter of the mode l, related to the grai n emi ssion efficiency, 
the same as used in IWvatt et all (|2012l ) and lLestrade et al.l 
(|2012t ). 

The spatial distribution of the dust is assumed to be 
disk-like, with a small opening angle. It is characterised by 
the dust's face-on optical depth (r), which is parametrised 
as a function of radius by one or more power laws. The dust 
can be arranged in a single or multiple belts, characterised 
by their radial location, width, position angle and inclination 
to the line of sight. The quality of each model is evaluated 
using a sum of squared image-m odel differences approach as 
described in [WyatteTal] (|2012l ). 
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Figure 7. The residuals, in units of significance of (image- model)/uncertainty, from the model fits with a single belt described in i|4.2l 
at 100 (left) and 160/im (right). The models provide a good fit to the observations, the residuals are at low levels (note the scale), with 
the black contours showing residuals of 3 — o". The white contours show the 3 — cr detection of the disc, as shown in Fig. [3] and Fig. [4] 




20 10 -10 -20 20 10 -10 

a cos 8 offset (") a cos 8 offset (") 



Figure 8. The same as Fig. but this time showing the residuals from the two belt fit. 

These models are defined by a large number of param- features. Differences may also exist between different sets 

eters, which means that the best fit model is not deter- of photometry due to calibration offsets. As can be seen in 

mined by undertaking a grid search of all possible param- the SED (Fig. |6} the spectrum suggested by the four far-IR 

eter spaces, rather by a combination of by-eye fitting and points is not well represented by the modified black-body, 

least-squares minimisation. While we do not claim that the thus, it may be suffering from these problems. While these 

models are unique, we do show that they are plausible lay- differences do not pose a major problem for SED modelling, 

outs of the debris disk structure. As we show below, we can they complicate image modelling since a small absolute off- 

reproduce the observed disk structure with several different set caused by an attempt to achieve the best fit for all data 

models, and therefore consider the uncertainties on Individ- can make an otherwise satisfactory image model to appear 

ual model-specific parameters relatively unimportant. very poor when compared with the ob s erved image. Our 

approach, therefore, follows [Lohne et al.l (|2012l ). who intro- 

Given that our modelling approach does not include juce small offsets at some wavelengths, so that the image 
detailed grain properties, it is likely that the true disk emis- models focus on fitting the radial distribution of the emis- 
sion spectrum is more complex than our simple modified ^^^^^ rather than the overall disc flux, which can be biased 
black-body, for example including relatively narrow spectral 
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due to the preceding factors. In this manner we could apply 
small modifications to each of the four far-IR points, how- 
ever, for simplicity we found that it was only necessary to 
apply such a modification to the PACS lOO^m point, using 
a factor Cioo- 

Considering firstly a single dusty ring, our best fit model 
is a wide belt extending from 20 to 220AU, with an opti- 
cal depth of r = 2.5 x lO^^r*^'^, a temperature profile of 
T = ^%7r-°'' I.e. fr = 1.1, Ao = 70Aim, /3 = 0.6 and 
Cioo ~ 0.95. These parameters are not particularly well con- 
strained, for example, from the SED fitting we already con- 
cluded that Ao and /3 are unconstrained, and reasonable fits 
to the data can be made for different values of the disc in- 
ner radius. However, a clear conclusion that a single, narrow 
belt does not fit this data is made. Better constrained are 
the inclination and position angle, for which we determine 
values of 59° from face-on and 142° , respectively. The values 
quoted here are similar to those found by fitting an ellipse to 
the emission in i]3.1l and we consider these to be constrained 
to within ±10°. A reduced of 0. 8 is calculated for t his 
model, in the manner described in IWvatt et all (|2012l '). 
The residuals once this model fit is subtracted from the ob- 
servations (Fig. O show that this model is a good fit to the 
data. The positive increase in optical depth with radius is 
unusual for debris disc models, however, we were unable to 
find a negative slope that provided a good fit to the data. 
This positive increase in surface density could be a real fea- 
ture, related to the dynamics or stirring mechanism of the 
disc. Alternatively, it could provide an indication that this 
structure is an incorrect interpretation of the real disc struc- 
ture. 

Given that such a wide single belt fits the data, the 
possibility that the system could contain multiple belts is 
worthy of investigation. We therefore relax the constraints 
and allow the model to include two dust belts. So that there 
is no increase in the number of parameters required in order 
to obtain a fit, we fix various parameters. We assume that 
both belts are lOAU wide, that the temperature dependence 
is T = 521r"°-^ (/t = 1), Ao = 70^im and that each ring 
has constant optical depth with radius. We find Cioo = 0.91. 
The best fit model has two belts, centred on 41 and 165AU, 
with optical depths of r = 2.7 x 10"'' and r = 1 x 10~^ and 
/3 = 0.7 and 1.0, respectively. The disc orientation, with a 
position angle of 145° and inclination of 60° from face-on, 
remains close to the original estimation. The residuals for 
this model are shown in Fig. \E\ and a value of = 0.7 
indicates that this is a good fit. 

Both of these models provide equally good fits to the 
data. They both have a common radial scale for the dusty 
material (from around 20AU to around 200AU) and an in- 
ner hole. The differences between these models reflects our 
lack of knowledge of the distribution of the emission within 
this region. In fact, the inner hole need not be completely 
empty and the observations are equally consistent with a 
stirring model in which the disc extends from near to the 
inner planet out to large radial distances, but is coUisionally 



^ Although the number of degrees of freedom (7 wei ghted values) 
used t o calculate this value may be misleading (see lWvatt et al] 
||2012| ')). a value close to 1 suggests a satisfactory fit. 



eroded from inside-out. In such a model the majority of the 
emission is still from the 20-200AU region. 

Planetesimals can be stirred as a result of the manner 
in which the protoplanetary g as disk was dispersed, the for- 
mation of Pluto-sized objects (|Kenvon fc Bromlev 20041'). or 



pertu rbations by a planetary companion (|Mustill fc WvattI 
120091 ). Given the known orbit for the inner planetary com- 
panion, we can estimate the radial distance out to which 
the planet could have stirred the disc. Using Eq.l5 of 
iMustill fc Wvat'3 (|2009l ). equating the stirring timescale, 
with the 2.5Gyr age of k CrB, the 2.1Mj planet could stir 
planetesimals out to around 70AU. The second, more dis- 
tant companion might inevitably stir planetesimals beyond 
this distance, depending on the properties of the orbit. In 
order to fit such a model to the observations, we consider 
that the surface density distribution of dust is split into two 
regions, an inner region where the density increases with 
radial distance to the star and an outer region in which 
it decreases. We then determine the radius of this change 
in surface density, which for these images occurs at 81AU, 
with the outer region extending to 290AU. Our best fit finds 
r = 5.8 X lO'V^-^ inside of 81AU and r = 6.6 x 10"^-^ 
in the outer regions, where the slope of the optical depth in 
the outer region was fixed at r~^, in both regions was fixed, 
with /t = 1.15. This model provides as good a fit to the 
data as the previous two models, as indicated by a x^ value 
of 0.8. Again, however, these parameters are not very well 
constrained. This model fits well with the planet stirring of 
planetesimals out to around 70AU by the known companion, 
but if our radial constraint is correct, the second companion 
must have orbital parameters close to the minimum possible 
(flpi = 7.3AU and Mpi = 0.5Mj) and orbit on a very low 
eccentricity orbit {Cpi < 0.06) in order that planet-stirring 
of planetesimals does not occur beyond ~ 80AU. 

To summarise, we have presented three plausible mod- 
els for the dusty material, as illustrated in Fig. (5] All three 
models reproduce the Herschel images and the SED. With- 
out further observations we have no means of distinguishing 
between these models, nor ruling out alternative models. 



5 DISCUSSION 

In this work we have presented new Herschel resolved im- 
ages of the debris disc around k CrB, alongside evidence 
for a second companion from Keck radial velocity data and 
upper limits on its mass from Keck AO im aging. Coupled 
with th e known planetary companion from I Johnson et al.l 
(|2008bl ). this allows us to constrain the structure of the k 
CrB planetary system. Although our knowledge of the k, CrB 
planetary system has grown significantly from these obser- 
vations, it is critically important to distinguish between the 
information that is well constrained from these observations 
and the more tentative conclusions that can be made. This 
will be discussed in the following section. 

By constraining the structure of the k. CrB planetary 
system, we have provided an example of a planetary system 
around an intermediate mass subgiant, or 'retired' A star, 
which in turn aids our understanding of the population of 
planetary systems around higher mass stars. In the second 
half of this discussion, we consider the impact of this study 
on our understanding of planetary systems in general. 
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Figure 9. A cartoon to illustrate the possible configurations of the k CrB planetary system. The top panel illustrates the possibility 
that the second companion lies interior to the observed dusty material, that lies either in a single wide belt or is split into two narrow 
belts, as described by the best fit models in ^4.21 The middle panel illustrates the possibility that the outer companion is in fact a binary 
and orbits between the two narrow dust belts. In the latter scenario the outer dusty belt would bo a circumbinary debris disc. The lower 
panel illustrates the stirring model, in which the rate of dust production is maximum at ~ 80AU. Diagram is not to scale. 



5.1 The Structure of the k CrB Planetary System 



Firstly, we consider the Herschel resolved images. One of the 
clearest conclusions that can be made regards the inclination 
of the circumstellar disc. The ellipsoidal nature of the excess 
emission can be clearly seen in the star-subtracted images 
shown in Fig. [S] The ellipse-like nature of the source sug- 
gests an azimuthally symmetric, circumstellar disc, viewed 
inclined to the line of sight. The inclination of this dust belt 
is reasonably well constrained from the images, at around 
60° from face-on, with a position angle of 145° . This in turn 
has implications for the second companion. If the planet and 
the dust disc formed out of the same proto-planetary disc 
and there have been no further interactions, it seems reason- 
able to assume that the inclination determined for the disc 
is also the inclination for the planet. This would mean that 
the planet's mass is a factor of ^.^g^, ~ 1.15 higher than 
the minimum mass of msini = 2.1Mj, with similar impli- 
cations for the mass of the unconstrained companion. The 
dotted line on Fig.[5]shows this increase in the planet's mass. 
It should, however, be noted that the disc and the known 
planet are sufficiently well separated that it is possible that 
post formation processes have altered their inclinations from 
coplanar. Although, during the the 2.5Gyr lifetime of k CrB 
, the secular perturbations from the known 2.1Mj compan- 
ion could have aligned the disc out to ~ 120AU, as suggested 



in (e.g. lAugereau et al.|[l999l : iKennedv et al.ll2012bl 'FI. This 
means that if the disc and planet were misaligned early in 
their evolution, the disc might appear warped at around 
~ 120AU, and our observations do not have sufficient reso- 
lution to detect this warp. 

The radial location and structure of the debris disc is 
not as well constrained as its inclination. In fact the main 
conclusion that can be made from these Herschel observa- 
tions is that there is clearly a dense population of dusty 
bodies orbiting k CrB. We suggest three structures that 
could explain the observations, although, based on the cur- 
rent data, we are unable to rule out other possibilities. Our 
first model is a continuous dust belt, extending from from 
20 to 220AU. The second is two distinct narrow dust belts, 
centered on 41 and 165AU. The third is a stirring model, in 
which the disc is coUisionally eroded from inside outwards. 
In this case, the dust production peaks at around 80AU. 
Currently our only source of further information regarding 
the structure of the planetary system comes from a consid- 
eration of the effect of the second companion on the dusty 
material. 

Our constraints on the orbit of the second companion, 
shown in Fig. O are not tight. We can, however, examine 
how it could fit into the three scenarios proposed. These are 
illustrated in Fig. (9] In all three cases, the second companion 
could orbit close to the inner companion, interior to the 

^ This is the distance at which particles would have processed 
once about the planet's orbital plane. 
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dusty material, be it in a single belt, two distinct belts, or 
distributed in a wide belt that is stirred by the planet(s). 
The companion could even be responsible for sculpting the 
inner edge of the inner belt. 

On the other hand, if the dusty material were to be split 
into two narrow belts, we consider the possibility that the 
second companion could orbit between the belts. In this case, 
the companion could be responsible for clearing the region 
of material. The outer belt at 165AU is sufficiently far from 
the second companion (that must lie within 70AU) that it is 
unlikely to be affected by it. The stability of a belt between 
36 and 46AU, however, is likely to be strongly influenced by 
the second companion. We note here, that our stated values 
for the exact radial position of the inner belt are not strongly 
constrained and that reasonable fits to the data can be made 
using different radial distributions of dust within the same 
region. Assuming that the companion is on a circular orbit, 
we can use the overlap of mean-motion reson ances, outlined 
by the chaotic zone {Sachaos = 1.3(^)2/^) (|Wisdomlll980l ') 
to estimate the size of the region in which particles would 
be unstable and thus, we anticipate to be clear of dust. The 
derivation for the size of this zone is only strictly valid for 
companion masses that are significantly less than the stel- 
lar mass, which is not the case for companions >> lOMj, 
nonetheless it remains true that higher mass planets must 
be further from the disc, if they are not to destabilise it. 
Using this limit, we find that if the companion were to be 
at the maximum radial distance from the star of 70AU, and 
therefore have a mass between ~ 42Afj and ~ 47Mj, a belt 
of dust could only be stable interior to this companion, if its 
outer edge were to be interior to 34AU. The values predicted 
for the inner belt in our model rule out such a scenario, how- 
ever, given the uncertain nature of our radial constraints on 
the dusty belts, the only conclusion that we can make is that 
if the second companion is to orbit between the two belts 
of dust, it must be on an orbit with semi-major axis close 
to 70AU. The limits on the mass of the second companion 
are such that in order for it to orbit between the two belts 
it must be sub-stellar, rather than planetary in natur^ll- lu 
this case, hi CrB would have an outer circumbinary debris 
disc, orbiting a main-sequence and brown dwarf binary, in 
which the main-sequence star has its own debris disc. 

To summarise, n CrB is orbited by significant levels of 
dusty material and at least two companions, at least one of 
which is planetary. The inclination of the planetary system is 
reasonably constrained with an inclination of 60° from face- 
on and a position angle of 145°. We present three plausible 
models for the distribution of dust, a single wide belt, two 
narrow belts or a stirring model. Further resolved imaging 
of this source in the infrared or sub-mm, combined with 
further radial velocity monitoring of hi CrB to constrain any 
outer planetary companions are required in order to be able 
to fully constrain the structure of this planetary system. 

5.2 K CrB in context 

The hi CrB planetary system is special, firstly, because it is a 
rare example of a system where both planets and planetesi- 

^ Using M > 13A/j, the deuterium burning limit, to define tlie 
maximum mass of a planet 



mal belts have been detected orbiting an intermediate mass 
star (m > I.4M0) and, secondly, because it is a unique 
example of a debris disc around a subgiant. hi CrB is the 
only > I.4M0 star with a giant planet inside of 8AU and a 
resolved image of a debris disc. 

The evolution of the star should not, on the subgiant 
branch, have had an effect on the planetary system, partic- 
ularly not the outer regions of the planetary region studied 
in this work. The main difference between this and most de- 
bris discs observed around main-sequence stars, is that the 
debris disc must have survived the en tire main-sequence life - 
time. K CrB has an age of 2.5 Gyr (| Johnson et al.ll2008al '). 
most of which was spent on the main-sequence. Models for 
the coUisional evolution of debris discs show a decrease in 
their fractional lumin osity with time, as c ollisions erode the 
material in the disc IWvatt et al.ll2007al '). An extension of 
these models, that also include the evolution of the star 
(I Bonsor fc Wvatt|[2O10l ). shows that this coUisional evolu- 
tion is the main factor that diminishes the detectability of 
debris discs around subgiants. A population survey for de- 
bris discs around subgiants is required to confirm whether 
K CrB is unusual in retaining its debris disc until the sub- 
giant branch. The fact that we detect a debris disc at this 
epoch implies that hi CrB did not suffer an event similar to 
the Late Ifeavy Bo mbardment that cleared our Solar System 
jBooth et al.ll2009l '). 

The main advantage of the subgiant nature of this star 
is that it enables radial velocity techniques to find plan- 
etary companions in a manner that would not have been 
possible whilst this star was on the main-sequence. This 
raises the question as to whether other main-sequence A 
stars that have debris discs similar to that of hi CrB may 
also have undetected planetary companions. For example, if 
hi CrB has a single wide belt, then it is not dissimilar to the 
single wide planetesimal belt of V ega (|Aumann et al.l [19841 : 
ISu et al.ll2005l : iMiiller et al.ll2010l '). which raises the question 
of whether Vega potentially has inner planetary companions, 
that cannot currently be detected. If on the other hand, the 
dust in K CrB is split int o two narrow belts, similarities can 
be se en with HR 8799 (|Su et al.1 12009| : iMarois et all |2008| . 

in that the dust belts are separated by planetary com- 
panion(s). In turn, this comparison could suggest that k, CrB 
may have further, undetected companions. 



6 CONCLUSIONS 

We have presented observational evidence for the structure 
of the K CrB planetary system, a unique example of a debris 
disc around a subgiant and a rare example of an interme- 
diate mass star, where both planets and planetesimal belts 
have been detected. Herschel observations show high levels 
of excess emission, evidence generally taken for the presence 
of a debris disc. Following detailed modelling of the Her- 
schel observations we suggest three possible structures for 
the dusty material. Either a single wide belt, extending be- 
tween ^ 20 and ~ 220AU, two narrow belts, centered on 
~ 40 and ~ 165AU, or a stirred belt, most probably stirred 
by the planetary companion(s), in which the dust produc- 
tion rates peak at ~ 80AU. Our best constraint is on the 
inclination and position angle of the disc, which we place at 
60° from face-on and 145° , respectively, to within ~ 10° . 
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Alongside the Herschel observations, we have presented 
evidence for the presence of a second companion to k CrB 
found in the continued radial velocity monitoring of this star. 
This is in a ddition to the m sin i = 2.1 Mj planet, at 2.8AU, 
detected bv I Johnson et al] (|2008al ). An upper limit on the 
mass of this companion was found by its non-detection in AO 
imaging taken with Keck. The details are shown in Fig. [5] 
In terms of our suggested models for the dusty material or- 
biting ft CrB , one possibility is that both companions lie 
interior to all the dusty material, potentially sculpting the 
inner edge of the inner belt. Alternatively, if the dusty mate- 
rial is split into two belts, we cannot rule out the possibility 
that the second companion could lie between these belts, 
giving K CrB an intriguing structure with both a circumbi- 
nary and circumprimary debris disc. 

As the first example of a planetary system orbiting a 
subgiant, a more detailed population study is required to 
determine whether or not n CrB is unusual, nonetheless, 
this work suggests that k CrB did not suffer any dynami- 
cal instability that cleared out its planetary system, similar 
to the Late Heavy Bombardment. As the first example of a 
> 1.4Mq star, with a giant planet interior to 8AU, where 
there is also resolved imaging of a debris disc, n CrB pro- 
vides a good example system from which to further our un- 
derstanding of planetary systems around intermediate mass 
stars. 
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